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to the polarization echo in glasses 
(extended) 

A. Borisenko 1 ' 2 '^ and G. Jug 2,3 

1 NSC Kharkiv Institute of Physics and Technology, Akademichna 1, 61108 Kharkiv, Ukraine 
2 Dipartimento di Fisica e Matematica, Universitd degli Studi dell'Insubria, via Valleggio 11, 22100 Como, Italy 
3 INFN-Sezione di Pavia and CNR-IPCF Unitd di Roma, Italy 
(Dated: January 19, 2013) 

Startling magnetic effects on the spontaneous polarization echo in some silicate glasses at low 
and ultra-low temperatures have been reported in the last decade or so. Though some progress in 
search of an explanation has been made by considering the nuclear quadrupole dephasing of tunneling 
particles, here we show that the effect of a magnetic field can be understood quantitatively by means 
of a special tunnel mechanism associated with paramagnetic impurities. For the Fe 3+ -, Cr 3+ - and 
Nd 3+ -contaminated glasses we provide reasonable fits to the published data as a function of applied 
magnetic field and temperature. 

PACS numbers: 77.22. Ch, 74.55. +v, 61.43.Fs, 71.38.Ht 



5h 



i 

C 

O 

o 



> 



O 



X 



Anomalous low-temperature magnetic properties of 
some silicate glasses have attracted a considerable ex- 
perimental and theoretical interest since the unexpected 
findings by Strehlow et al. jjj . Since the low-temperature 
behavior of glasses is believed to be governed by tunnel- 
ing systems (TS) (see e.g. 0), a lot of theoretical effort 
has been applied to account for the coupling of TS to a 
magnetic field. All the theoretical models that have been 
proposed so far may be grouped into either " orbital" 
or "spin" type according to the assumed mechanism 
for the TS coupling to a magnetic field. 

A model of paramagnetic tunneling systems (PTS), be- 
longing to the "spin" group, has been recently proposed, 
in collaboration, by one of the present authors [7|. It 
considers an effective charged particle with an intrinsic 
spin S, tunneling between the coordinate states (wells) \l) 
and |r) of a ID double- well potential. In each coordinate 
state the spin is quantized along the axis z and the only 
spin projections allowed are S z = ±S . Therefore, one 
can treat both coordinate and spin states of such PTS 
in the spin- 1/2 formalism. Let two sets of Pauli matri- 
ces <r and t together with the corresponding 2x2 unity 
matrices I CT and 1 T refer to the PTS coordinate and spin 
states, respectively. Then the Hamiltonian of such PTS 
is as follows [a|: 



The eigenvalues of the Hamiltonian (|TJ) read: 
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where h is the energy difference (asymmetry) between 
the states |Z) and |r), u is the energy difference between 
the states \±S) in each well and the tunneling matrix 
elements are 5 = A cos (a/2), d = A sin (a/2), A being a 
tunneling matrix element between the states \l) and |r) 
of an ordinary (spinless) TS, and a an angle between the 
quantization axes zin and . 



Ux,A = t\\J(\u\ + VT^~P) 2 + d\ 

U2,z = t\ \j(\u\- y/tfTS*) 2 + <P. (2) 

Previously one of the present authors has shown that 
under certain temperature and field conditions the PTS 
model is applicable to a special type of paramagnetic 
impurity- hole complexes [FeO^] , presumably present in 
the Fe 3+ -contaminated silicate glasses 0. Under cer- 
tain model approximations, this resulted into a good 
agreement with experimental data for the specific heat 
and dielectric susceptibility of several types of iron- 
contaminated glasses as functions of temperature and ap- 
plied magnetic field [8(. 

The idea of existence of the [FeO^] complexes in sil- 
icate glasses is fully compatible with the explanation, 
given by Castner et al. Q for the broad EPR absorption 
line at a w 6 , detected already in the pioneering work by 
Sands 10]. It is also supported by the evidence for exis- 
tence of the, to them closely related, [FeOJ^ complexes, 
coming from interpretation of EPR spectra in quartz [ill ] . 

In the present paper we consider the paramagnetic 
impurity- hole complexes of the general type [XO4] , 
where X is a 3-valent substitution paramagnetic impurity. 
A schematic structure of the [XO4] complex is shown in 
FigQ] The ion X 3+ substitutes Si 4+ in the silica ele- 
mentary cell [SiO^] . Then one of the neighboring oxygen 
ions traps a hole to become O - (meaning the 2_ +hole 
complex) and to compensate the local charge defect [HI ■ 
This structure is isomor phic to the well-known [AIO4] 
color center in a-quartz [l3| . In the tetrahedral cage of 
oxygen ions, with the hole localized at one of them, the 
ion X 3+ is subject to a crystal field with an approxi- 
mate 3-fold axial symmetry along z, this being the local 
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FIG. 1. (Color online) 2D schematic structure of the 
impurity- hole complex [XO4] . 



basis axis, directed from X 3+ towards O - . A gradient 
of this crystal field couples to the electronic quadrupole 
moment of X 3+ and causes the quadrupole splitting of 
its spin states. The spin J of the paramagnetic ion also 
couples to the spin j = 1/2 of the hole and to the ex- 
ternal magnetic field B. Other interactions, such as the 
nuclear quadrupole, dipole electron-nuclear and nuclear 
Zeeman ones, are much smaller and may be neglected. 
Then, neglecting a possible nonaxiality of the local crys- 
tal field and the dipole interaction, the spin Hamiltonian 
of the [XO4] complex in either coordinate state may be 
considered in the following form: 
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5oj + 9x J • B + D z j z J z + Q z Jl . (3) 



Here /3 is the Bohr magneton, go « 2 and gx are the 
Lande factors for the ions O - and X 3+ , respectively, and 
D z < and Q z < are the dipole and quadrupole inter- 
action constants, respectively. 

In the experimentally relevant range of magnetic fields 
B ^ J\D z \/go/3cos(a/2) ~ 0.1...1 T, the ground state 
of the Hamiltonian (J3j) is a quadruplet of spin states 
|±J, ±1/2). The first (non- vanishing) order corrections 
from the Zeeman and dipole interactions to the ground- 
state energy level Q Z J 2 are as follows: 
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In the temperature range T < \D z \/k-Q ~ 0.1 K the 
doublet of levels 1 and 2 is the only thermodynamically 
relevant one. In this case the Zeeman splitting value is 
3: 



u^2f3(± go + Jgx)B cos (a/2). 



In the temperature range \D z \/k^ < T < \Q z \/kB ~ 
10 K the whole quadruplet of levels (HJ is thermodynam- 
ically relevant. Strictly speaking, the PTS model in its 
simplest form ([TJ is inapplicable to the case of more than 
two relevant spin states. Nevertheless, one can hope that 
it gives a qualitatively adequate physical picture even in 
this case. At B = the only gap in the spectrum §4§ 
is due to the dipole interaction: u (0) = J\D Z \. For the 
purpose of illustration of the model's results at finite val- 
ues of B, as an heuristic assumption, we take the gap 
between the outside levels 1 and 4 as the Zeeman split- 
ting value 
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Now, if one considers tunnel transitions between the 
states |Z) and |r) in FiglTJ the PTS model becomes appli- 
cable to this physical object. 

In this paper we study the PTS contribution to the am- 
plitude of the two-pulse (spontaneous) polarization echo 
and compare the results to experimental data for several 
glasses with Fe 3+ , Cr 3+ and Nd 3+ paramagnetic impuri- 
ties. 

The essence of the spontaneous polarization echo phe- 
nomenon in glasses is as follows. At t = the first mi- 
crowave pulse with frequency ujq and duration T\ creates 
a macroscopic coherent superposition of a subsystem of 
TS with energy gaps At/ ~ huQ. This macroscopic co- 
herence quickly vanishes due to the wide distribution of 
the local TS parameters, but the microscopic (individual 
phase) coherence remains much longer. At t — t w — T2/2 
the second pulse with duration ti reverses the temporal 
evolution of the TS. Then at t w 2t w the macroscopic co- 
herence is spontaneously restored in the form of an echo. 

An external harmonic electric field with an amplitude 
Fo couples to the PTS dipole momentum operator P = 
7;Pol T <8> <J Z and adds a perturbation term V (t) = Fo • 
p e -^ t + ( Fq . p)+ e ^ t into the PTS Hamiltonian ©. 
In this case one can consider the following form of the 
PTS wave function: 



PTS 



(t) = a n (t) exp (-iU n t/h) \ip n ) 



n=l 



(7) 



where \ip n ) are the eigenstates, corresponding to the 
eigenvalues U n © of the PTS Hamiltonian JT]). The 
coefficients a n (t) obey the normalization condition 

E \a n (t)\ 2 = l. 

n=l 

If only one gap in the PTS energy spectrum ((2]) is close 
to resonance: 

Wmfc = {U m - Uk) /h = UJQ + e m k, \£mk\ < ^0, (8) 

the corresponding coefficients a m and ak are governed by 
the following equations (ljj : 



(5) 



ida m /dt = r\ mk exp (is m kt) ak 
ida k /dt = rf mk exp (-ie mk t) a ri 



(9) 
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where r] mk = (tp m \ Fo • P \<Pk) /fi is nonzero only in pres- 
ence of the microwave field. 

For the purpose of solution of our problem it is conve- 
nient to change to new variables: 

b» (t) = a n (t) exp {-iU n t/K) . (10) 



In these notations the system of Eqs. ([9]) reads: 

idb m /dt = rimk exp (ie mk t) b k + U m b m /h 
idb k /dt = i]* nk exp (-ie mk t) b m + U k b k /h. 

At t = Ti (immediately after the first pulse) the solution 
of Eqs. (jlip is as follows: 



n mk b k (0) cos(q nfc ri) +isin(n mfc TL)[(emfc/2)bfc(0) ~ < lfc ^(0)]e~ t£ "' fcri/2 

£l mk exp[iU k Ti/H] 

^mfcfem(0)cos(» mfc ri) - ism(n mk T 1 )[(e mk /2)b m (0) + y m kb k (0)}e l£ ™ kTl/2 

exp[iU m Ti/h] 



Here fl mk = y\r] mk \ 2 + £„ fc /4 is a generalized Rabi fre- 
quency. 

During the first period of free evolution t\ < t < t w — 
t 2 /2 the non-diagonal coefficients rj mk are zero, therefore 
just before the second pulse the solution of Eqs. (fTT]) 
reads: 

b k (t w - r 2 /2) = b k (n) eiq)[-iUh(t w - r 2 /2 - n)/H], 

b m (t w - T2/Z) = b m (rx) exp[-iU m (t w - r 2 /2 - n)/h]. 

(13) 

Bearing in mind, that only those PTS make a contribu- 
tion to the echo signal, which undergo a transition (emit 
or absorb a photon) during the second pulse, one can 
readily find the solution of Eqs. (JTTJ) at t > t w + t 2 /2 (in 
the period of free evolution after the second pulse): 

bh(t -t w - r 2 /2) = 

-ir,* mk sm(n mk T 2 )b m (t w - T 2 /2)e- le ™ kT ^ 2 

Q. m k exp[iU k (t — t w + r 2 /2)/h} 
b m (t -t w - r 2 /2) = 

-iT)mk sm(n mk T 2 )b k (t w - T 2 /2)e l£mkT2/2 



exp[iU m (t — t w + t 2 /2)/H] 



(14) 



The solution of Eqs. (JTTJ) in the form Eqs. ([12]) - ([13]), 
with the initial conditions 



M0)| 2 = k(0)| 2 =nr(T) 
4 

/ (■■; \ / < 

exp 




(15) 



given by the Boltzmann population of the PTS levels 
([2]) before the first pulse, can be used to find the time 



dependence of the PTS polarization: 

^PTS (t) = (1-PTS (*)|P0 ' P I^PTS (t)> /Fo- (16) 

Neglecting all possible dynamic-dephasing mechanisms 
such as the interaction with phonons and with neighbor- 
ing TS and omitting the terms odd in rj mk (they van- 
ish due to the isotropic angular distribution of the PTS 
dipoles po), one can derive from Eq. (|16p an explicit 
formula for the echo amplitude (envelope): 

A PTS (0 =ioE \r)mk\ (nT ~ <?) H mk (f ) , (17) 



where 



H mk (*') = sin 2 (O m fcT 2 ) sin (2fi TOfe Ti) cos e mk t' + 
(£ mfc /il mfe )sin 2 (fi mfc ri) sin£ mfe i' |?7mfc| 3 /^;L- (18) 

Here ^4o cx fi/-Fb is a prefactor and t' = t — 2t w +t±. In 
Eq. ([TS]) we have omitted the terms odd in e, since they 
vanish after integration over this variable (see Eq. ([20]) 
below) . 

An average echo amplitude is a result of integrating 
Eq. (jTTJ) with a distribution function for the double- 
well potential parameters f (h, A, a). In practice, it is 
convenient to change one of the integration variables to 
e, e.g. h — > e. As an approximation, we integrate only 
the "fast" function of e Eq. (|18p. The remaining ones 
are taken at the resonance hypersurface: 



lk (h,A,a,u) = 0. 



(19) 



This is valid since Eq. (1181) is nonzero in the narrow 
region \e mk \ < \i] mk \ oc |p • Pq| /H < uiq. In this way 
one is able to represent the average echo amplitude as 
follows: 
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f (h, A, a) / de mk H mk (t') 



(20) 



where h mk (A, a, w) are the real roots of Eq. (fl9)). 

An integrated echo amplitude is a result of integrating 
Eq. HI over i': 



dt'Ap T s (*') . 



(21) 



To compare the results of the PTS model to the available 
experimental data we use the following model distribu- 
tion function for the double-well potential parameters: 



f(h,A,a) 



1 



2ir6h 



exp 
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2nSa 



exp 



(a - a ) 2 
26a 2 



(22) 



with 5h/k B = 230 K, A /k B = 40 mK, A max /fc B = 4K 
and the parameters ao and da depending on the sub- 
stitution impurity as listed in Table |U Eq. (|2"2")l is an 
analytic function of A: it is zero at A = 0, then passes 
through a maximum at A = A and falls asymptotically 

as 1 / A at A ^ Ao , similar to the standard TS distribu- 

r~ii 

tion function (see e.g. |2|). 

As it was stated above, the temporal attenuation of 
the echo signal is out of the scope of our considera- 
tion. Ludwig et al. fisj ] found experimentally that in the 
barium-aluminosilicate glass (BAS), studied at 1 GHz 
and 10 mK, the integrated echo amplitude depends on 
the waiting time exponentially: A lnii cx exp (— t w /r), with 
the relaxation time r being practically independent on 
the value of applied magnetic field. Therefore, one can 
conclude that the relative value of the echo amplitude, 
e.g. A{B) /A(0), should be independent (roughly, at 
least) on the experimental value of the waiting time t w . 
In this way one also gets rid of the main temperature 
dependence of the echo signal due to the temperature 
variation of the relaxation time r, with the only residual 
temperature dependence coming from the initial equilib- 
rium levels' population n° q (T), given by Eq. p3|) . 

From the numerical studies of Eq. (|20|) as a function of 
applied magnetic field at t' = we find that it reaches a 
minimum at some value -B m i n and that within our accu- 
racy Apxs (-Bmin) ~ 0- I n terms of the Zeeman splitting 
u the minimum occurs at zt m ; n ~ hu Q . 

With the above speculations in mind, in 
Fig. [2] we re-plot experimental data 1_7[ for 
the integrated echo amplitude in the Fe- 
contaminated glasses Duran and BAS in the form 
[A- mt (B) - A int (B min )} I [A int (0) - A int (B min )] (to 
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FIG. 2. (Color online) Integrated echo amplitude vs. B. The 
points stand for experimental data [13] • The curves are calcu- 
lated from the present theory. The data are shifted arbitrarily 
along the vertical axis. Details as given in the text. 



subtract the magnetic-independent contribution, coming 
from the ordinary TS) and compare them to the function 
A int (B) /A- mt (0), calculated from Eq. ([21]), with the 
Zeeman splitting calculated from Eq. ([5]) and with the 
parameters, corresponding to the [FeO.4] complex, as 
listed in Table U 

In Fig. [3] we re-plot experimental data 
echo intensity in the Cr- and Nd-contaminated glasses 



18( for the 



in the form 1(B)/ 1(0) together with the function 
^pts ( B ) Mpts (°)> calculated from Eq. ® at t' = 0, 
with the Zeeman splitting calculated from Eq. ([5]) and 
with the parameters for the corresponding complexes as 
listed in Table |U In the case of the Nd-contaminated 
glass, for the purpose of fitting, we added a constant to 
Eq. ff20|) . which accounts for the magnetic- independent 
contribution, presumably coming from the non-magnetic 
TS, probably of the [AIO4] or [BO4] types. 

Yet another interesting phenomenon, reported [Til ] for 
the Cr-contaminated glass, is a dynamic enhancement 
of the echo intensity during a rapid decrease of the 
magnetic field from a constant value in the vicinity of 
B m ax ~ 0.3 T, corresponding to the maximum of the 
echo intensity (see FigEJ). In the framework of the 
present theory we provide the next qualitative explana- 
tion for this phenomenon. Since the decrease of the mag- 
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TABLE I. Fitting parameters for the different [XO4] complexes. 



Complex 


J 


gx 


ao 


6 a 


J\D z \/k B , (K) 


|po • F n/h 


Po ■ Fo \r 2 /h 


[FeOJ 


5/2 


2 


2.6 


0.1 




0.819 


1.638 


[Cr0 4 ]° 


3/2^ 


T 


2.07 





0.354 


2.235 


2.235 


[NdOJ 


3/21 


2 a 


2.42 


0.1 


0.22 


1.117 


1.117 



a These values correspond to the zero orbital momentum of X 3+ and, therefore, to the minimal ionic radius to minimize the positive size 
misfit of X 3 + compared to Si 4+ . 
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FIG. 3. (Color online) Intensity of the echo signal vs. B. 
The points stand for experimental data [18j. The curves are 
calculated from the present theory. Details as given in the 
text. 



TS is already damped by phonons. 

It seems reasonable that the echo response in the 
glasses free from para mag netic impurities, such as the 
borosilicate glass BK7 |l7| and deuterated glycerol [20j j . 
can be explained in the framework of the quadrupole 
model [6j. The limited applicability of the quadrupole 
model explains its success in reproducing the experimen- 
tal data only for BK7 @,[2l[ and glycerol [Hj]. 

It is tempting to assume, that TS of the same type, 
but originating from the non-magnetic substitution im- 
purities, such as the [AIO4] or [BO4] ones, also give 
their contribution to the low-temperature properties of 
multisilicate glasses and may be responsible for some de- 
viations of experimental data from the predictions of the 
standard tunneling model. 

In summary, we have demonstrated, that with the 
assumptions made, the idea of application of the PTS 
model to the [XO4] paramagnetic impurity- hole com- 
plexes results into a systematic reasonable agreement 
with the experimental data from the spontaneous po- 
larization echo in different silicate glasses with different 
paramagnetic impurities and in a wide range of frequen- 
cies and temperatures. 



nctic field leads to the decrease of the Zccman splitting u 
(see Eq. ©), the population of the PTS levels turns into 
a non-equilibrium one, namely, the excited levels become 
under-populated. The equilibration process goes through 
the absorption of phonons by the PTS. This "pumping" 
of the energy from the phonon subsystem into the PTS 
one leads to the decrease of the sample temperature T. 
This process is analogous to the well-known adiabatic de- 
magnetization cooling of paramagnetic samples. Then, 
the decrease of the sample temperature T leads to the 
slowing down of the phonon attenuation of the echo sig- 
nal and, therefore, to its enhancement. 

It should be mentioned that, from the studies of the 
dielectric relaxation phenomena in smoky quartz (l9| . 
containing the [AIO4] complexes, one can expect the 
characteristic relaxation rate due to the interaction with 
phonons for the [XO4] complexes to be much slower than 
for the ordinary TS. That is why in the echo experiments 
they manifest themselves even in the liquid helium tem- 
perature range, where the contribution from the ordinary 
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